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Abstract. We present three-component wind models for carbon rich pulsating AGB stars. In particular we study the effects 
of drift in models of long-period variables, meaning that the dust is allowed to move relative to the gas (drift models). In 
addition we investigate the importance of the degree of variability of the wind structures. The wind model contains separate 
conservation laws for each of the three components of gas, dust and the radiation field. We use two different representations for 
the gas opacity, resulting in models with different gas densities in the wind. The effects which we investigate here are important 
for the understanding of the wind mechanism and mass loss of AGB stars. This study is hereby a necessary step towards more 
reliable interpretations of observations. We find that the effects of drift generally are significant. They cannot be predicted from 
models calculated without drift. Moreover, the non-drift models showing the lowest mass loss rates, outflow velocities, and the 
smallest variability in the degree of condensation do not form drift model winds. The wind formation in drift models is, except 
for a few cases, generally less efficient and the mass loss consequently lower than in the corresponding non-drift models. The 
effects of drift are generally larger in the more realistic models using that representation of the gas opacity which results in 
lower densities. The outflow properties of these models are also - for all cases we have studied - sensitive to the period of the 
stellar pulsations. A check of the mass loss rates against a (recent) fit formula shows systematically lower values, in particular 
in the more realistic models with a low density. The fit is in its current form inapplicable to the new models presented here. 
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1. Introduction 

Stellar winds of asymptotic giant branch (AGB) stars are be- 
lieved to form as a consequence of an efficient dust conden- 
sation in the extended atmosphere of these stars. Dust grains 
absorb radiative momentum from the highly luminous star and 
are thereby pushed away from it. While moving away from the 
star, they collide with gas particles and the momentum is trans- 
ferred to the rest of the matter, forming the stellar wind. Mass 
loss and its history is not only important to the AGB star itself. 
It is also relevant for its neighborhood, the formation of the 
circumstellar envelope and in a longer perspective, the galactic 
chemical evolution (e.g. Willson 2000). 

Since Bowen (1988) published his pioneering paper on 
time-dependent wind formation in cool M-type giants consid- 
erable progress has been achieved with models describing dif- 
ferent physical aspects in more detail. Fleischer et al. (1992) 
created a C-rich wind model using a self-consistent descrip- 
tion of dust formation. Fleischer et al. (1995) and Hofner et al. 
(1995) studied models without pulsations and described the 
dust-induced (external) kappa mechanism. The influence of 
molecular opacities for the wind driving efficiency were stud- 
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ied by Hofner et al. (1998, henceforth HJLA98) and Helling 
et al. (2000). The models in this article differ from the ones pre- 
sented in the papers mentioned above as drift between gas and 
dust is taken into consideration. So far this phenomenon has al- 
most exclusively been studied in stationary winds (e.g. Kriiger 
et al. 1994 (henceforth KGS94), Kriiger & Sedlmayr 1997, or 
Liberatore et al. 2001) or non-pulsating stars (e.g. Simis et al. 
2001 or Sandin & Hofner 2003 (henceforth Paper I)). 

In the first article in this series. Paper I, we presented a 
method on how to describe drift in a C-rich wind model. The 
work is based on the model presented in Hofner et al. (1995) 
and both aspects of physical and numerical details were cov- 
ered therein. The main objective in Paper I was to study the 
basic effects of drift on the wind structure; by comparing the 
results with position coupled (PC; i.e. non-drift) models. In par- 
ticular we studied different approximations of the term describ- 
ing the collisional gas-dust momentum transfer, the drag force. 
Processes that could make the interpretation of the results more 
difficult were excluded. Consequently we did not describe stel- 
lar pulsations in the models of Paper I. And the positive effects 
of atmospheric levitation by pulsation were hereby excluded, 
requiring a high carbon/oxygen ratio and a luminous central 
star. 
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The main result of Paper I was that the dust in drift models 
tends to accumulate to the dense regions behind shocks, which 
form as a consequence of the dust-induced kappa-mechanism. 
This accumulation of dust does not occur in PC models. While 
the differences were obvious in the spatial data we could not 
find any certain changes - between drift and PC models - in 
quantities which are temporally averaged at the outer boundary, 
such as the average mass loss rate and the average terminal 
velocity. The reason why we considered the averaged results of 
the drift models uncertain was that the numerical accuracy of 
the momentum flux terms in the drift models never was better 
than first order in Paper I. 

The final goal is to replace the PC models with drift models, 
if it turns out to be necessary for a specific purpose. However, 
before we can decide on this matter it is necessary to find out 
what effect the drift has on the models, and whether there are 
clear correlations of the drift to other quantities. Is it possible 
to describe the effects on the wind structure without explicitly 
including the drift in the calculations? Such a treatment would 
be time-economically valuable since the larger variability in 
the dust component caused by drift inevitably slows down the 
calculations. 

With the work presented in this article we improve the 
model presented in Paper I further The numerical accuracy 
in the flux terms of the partial differential equations is now of 
second order for all models. We here include a description of 
stellar pulsations to study long-period variables. A commonly 
used representation of the gas opacity in AGB wind models is 
a constant number introduced by Bowen (1988). For compari- 
son with the corresponding models we calculate models using 
this representation too. However, since its use results in unre- 
alistic density structures (cf. Paper I, Sect. 4.1) we also calcu- 
late models based on molecular data (using Planck mean ab- 
sorption coefficients; HJLA98). The selection of sets of model 
parameters is based upon a sample where we expected that 
drift effects could be important; the models in the sample are 
most closely related to the models presented by Hofner & Dorfi 
(1997, henceforth HD97) and HJLA98. The only related study 
of time-dependent dust-driven winds - we are aware of - that 
self-consistently treats dust formation and drift is presented by 
Simis (2001, ch. 4). However, the models presented therein 
are calculated using assumptions that are different from ours 
in several respects (cf. Paper I, Sect. 4.2). 

The physical assumptions and the model characteristics of 
the current work, compared to Paper I, are specified in Sect. 2. 
The parameter-relations and the criteria we use in the selec- 
tion of the sets of model parameters are given in Sect. 3, which 
also covers the modeling method. The results of the study are 
presented and discussed in two separate sections. The output 
of the models in terms of average properties of the wind and 
their inter-relations, are given in Sect. 4. Four topics associated 
to the new models and the effects of drift are discussed in fur- 
ther detail in Sect. 5. The topics are: the influence of the piston 
period on the properties of the wind; the importance of a "tem- 
poral variability" for the formation of winds and for discerning 
winds of different character; a discussion on the comparison 
with a mass loss fit formula; and the influence of heating by 
drift on the wind structure. The conclusions are given in Sect. 6. 



Table 1. Glossary of used symbols and units 



sym. 


unit 


■ ■ 
description 


D 
r 


g cm""* 


gas density 


Pa 


g cm"^ 


dust density 




cm"^ 


gas number density 




cm"^ 


dust number density 


u 


cms"' 


gas velocity 


V 


cms"' 


dust velocity 


J 


ergcm"'^ s"' 


zeroth moment of the 






radiation field 




cm"^ 


moments of the grain size 






distribution; < j < 3 


Vd 


cms"' 


drift velocity 






speed ratio 






drag coefficient 


£ 




fraction of specular collisions 


CT 


cm^ 


gas-dust geom. cross section 


/dmg 


gcm"-s"^ 


drag force 


Adrift 


erg cm"' s"' 


heating by drift 




cm 


mean grain radius 


P 


dyn cm"^ 


gas pressure 


T'i 


K 


gas temperature 


Ti 


K 


dust temperature 


5g 


ergcm"^ s"' ster"' 


source function of the gas 




g"' cm^ 


(gray) gas opacity 


Kd 


-1 2 
g cm 


(gray) dust opacity 


c 


cm s"' 


speed of light 


kB 


ergK"' 


Boltzmann constant 


G 


dyn g"^ cm^ 


gravitational constant 


0"B 


ergcm"^ s"'K"'* 


Stefan-Boltzmann constant 


Quantities specified at the outer boundary 


M 


Moyr"' 


mass loss rate 


Moo 


cm s"' 


terminal velocity 


ycond 




degree of condensation 


Td 




dust optical depth 


cr, 




standard deviation 


Grain properties 




Ko = tii 






= {rd)/ro; ro [cm] — 


monomer radius 




= <A>/4;rrg;<A>[cm 


-] — mean grain surface area 


K3/K0 


= (A'); (N) — mean 


grain size 



2. Physical assumptions and model 
characteristics 

The stellar wind can physically be separated into three different 
interacting components. All but a small fraction of the matter 
is assumed to form the gas component. The remaining mat- 
ter, present in dust grains, forms the dust component. The third 
component is the radiation field. 

In our wind model each of the three components is de- 
scribed by coupled conservation equations covering radiation 
hydrodynamics, time-dependent dust formation and interac- 
tion between all components; we call it the RHD3-system 
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(Radiation HydroDynamics, Dust and Drift). The physical sys- 
tem and how it is solved is described thoroughly in Paper I. In 
short the equation of continuity, the equation of (internal) en- 
ergy and the equation of motion describe the gas component. 
The dust component is assumed to be made of spherical par- 
ticles in the form of amorphous carbon and is described with 
the dust equation of motion and four moment equations treat- 
ing the formation and destruction of dust grains; the moments 
(K0-K3) describe certain properties of the grain size distribu- 
tion function. Currently dust formation is treated considering 
the processes of nucleation, growth, evaporation and chemical 
sputtering. The radiation field is represented by two moment 
equations of the radiation transfer equation. The two moments 
are the frequency-integrated zeroth and first moments of the ra- 
diation intensity, representing the radiative energy density and 
the radiative energy flux, respectively. All symbols and quan- 
tities used in this article are (for a quick reference) given in 
Table 1. 

In the three following subsections we look at processes 
where the physics used in Paper I is improved further. 

2.1. Treatment of the gas-dust interaction 

In this article we primarily present and discuss the results of 
drift models - i.e. models where the dust is allowed to drift 
with respect to the gas - but for comparison we also present 
PC models that a priori assume complete momentum coupling 
(and a tight mechanical coupling of the dust to the gas, pre- 
venting drift). In drift models the drift velocity vq is defined 

by, 



I'D = V - M 



(1) 



where v is the dust velocity and u the gas velocity. All three 
processes of momentum, energy and mass transfer between the 
gas and dust phases are affected by drift. 

The description of the momentum transfer between the gas 
and the dust is contained in the drag force which is calculated 
taking care of both the geometrical and thermal properties of 
the dust and the thermal properties of the gas. 



/drag = o-(Ko,Ki)pna 



''dCd(vd, T„, Td, e) 



(2) 



Here cr is the gas-dust geometrical cross section; Kq and Ki 
are the first two moments of the grain size distribution func- 
tion; p is the gas density; is the dust number density (i.e. 
Ko); Cd flie drag coefficient; Tg the gas kinetic temperature; Td 
the dust temperature; and s the fraction of specular gas-dust 
collisions (0 < e < 1). The drag coeflicient we adopt is C^^ 
given in Eq. (23) in Paper I. The exact form of the gas-dust 
geometrical cross section, the dust properties, and the assump- 
tions used for the dust temperature are given in Paper I. The 
term in Cd that accounts for diffusive collisions becomes sig- 
nificant at low drift velocities and large dust/gas temperature 
ratios, see Eq. (10) in Paper I. However, in Fig. 4c (in Sect. 4.3 
Paper I) it is seen that the same term never accounts for more 
than 40% of Cd and the calculated wind structure is more or 
less independent of the type of collisions. Consequently, with 



physical conditions similar to those described in Paper I we as- 
sume fully specular collisions (e = 1). 

The transfer of internal energy (dependent on drift) be- 
tween gas and dust can be split in two terms. On the one hand 
heat is transferred to gas particles that are accommodated on 
the dust grain surface in diffusive collisions. For a complete 
derivation of this term see e.g. Schaaf (1963, and the equa- 
tions leading to Eq. (7)); KGS94 (Eq. (14)) discuss an approx- 
imative form. On the other hand energy of the bulk motion is 
converted into internal (thermal) energy as gas particles which 
preferentially come from one direction (as in the case of drift) 
are reflected in random directions when hitting a dust particle. 
Calculations involving the details of a two-body collision show 
that this heating (by drift) is given by. 



1 +e 



Adrift 



-I'D/drag 



(3) 



(e.g. KGS94 and Draine 1980, 1986). KGS94 found heating 
by drift to be a significant factor in the energy balance of the 
gas in stationary drift models with a simplified description of 
the radiation field and the dust. In our time-dependent models 
we find that the impact of this heating on the wind structure is 
low; cf. Sect. 5.4. In contrast to Paper I, where g'drift was not 
included in any calculations, it is in this article included in all 
calculations presented in Sects. 4- &5.1. 

The mass transfer (represented by dust formation and de- 
struction) is also affected by drift, partly through modified grain 
growth efficiencies, partly in non-thermal sputtering that oc- 
curs frequently if the drift velocity exceeds about 40 km s ' (cf. 
Kriiger & Sedlmayr 1997, and references therein). In the re- 
sults of Paper I we found that a significant fraction of the dust 
in drift models accumulates to the regions behind gas shocks, 
the drift velocity was also found to stay below 30kms"'. The 
latest calculations presented in this article do not significantly 
alter the validity of the conditions behind these assumptions. 
Currently we neither include sputtering due to drift nor mod- 
ifications of the dust formation rates to account for drift. The 
plan is, however, to include them in the calculations and study 
their influence on the wind structure closely in the near future. 

2.2. Treatment of ttie gas-radiation field interaction 

Starting with Bowen (1988) it has been customary to assign the 
gas opacity a (low) constant number as a primary estimate. 



(4) 



It was done on the basis that reliable opacities at that time 
were not available and the value given here was thought to be 
representative in the region of the wind where the gas opac- 
ity was believed to play a role. To mention just a few, Bowen 
(1988), Fleischer et al. (1992), Fleischer et al. (1995), Hofner 
et al. (1995), HD97, Winters et al. (1997), Winters et al. (2000, 
henceforth WOO), and Paper I all present time-dependent dust- 
driven wind models that adopt this opacity. 

The constant value on the gas absorption coefficient (Kg) 
is independent of both the thermodynamical conditions and the 
chemical composition and its use results in unrealistically large 
densities in the photospheric layer It was replaced with Planck 
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mean absorption coefficients based on the SCAN molecular 
data (J0rgensen 1997) by HJLA98 and Helling et al. (2000). 
This improved description of the gas opacity avoids the extra 
"free parameter" (given in Eq. (4)), is a function of the chemi- 
cal composition, and allows for the calculation of a more real- 
istic density structure. 

The difference between structures of models using either of 
the two described opacities is in general significant, and the re- 
sulting winds show very different behaviors. In this article we 
both calculate models using a constant gas opacity (constant- 
opacity models) - to compare with the results of previous mod- 
els in the literature - and models using an opacity based on 
Planck mean absorption coefficients (Planck mean models). 

2.3. Modeling stellar pulsations 

To describe the effects of stellar pulsations on the atmosphere 
we use a sinusoidal radially varying inner boundary, located at 
about 0.91 /?« (above the region where the /c-mechanism sup- 
posedly originates). Since the radiative flux is assumed to be 
constant through the innermost mass shell the radial variation 
implies a varying luminosity as well. The above description of 
the pulsations is not self-consistent in the sense of a pulsation 
model, but it provides the wind with the dynamical effects of 
pulsations, i.e., a levitated atmosphere and strong shocks (cf. 
HJLA98 and references therein). The period of the piston is in 
the previous work of e.g. HD97 and HJLA98 specified through 
a period-luminosity (P-L) relation derived and fitted for M-type 
Miras of periods P shorter than 420 days by Feast et al. (1989), 

log P = 2.5/3- log (L/Lo)- 1.84/3 (5) 

where P is given in days. The piston amplitude is a "free pa- 
rameter", see Sect 3. 

The influence of different pulsation periods on the wind 
structure is studied by e.g. WOO and HD97. The former authors 
concluded that a long period allows for a stronger levitation of 
the atmosphere and more favorable conditions for dust forma- 
tion. However, these effects are counteracted by a more efficient 
radiation pressure on the dust that dominates the dynamics and 
cancels the improved efficiency. Similarly HD97 found the ef- 
fects of altered periods to be insignificant in two of their PC 
models. 

We want to see if this conclusion is valid for drift models 
and therefore complement the models calculated using Eq. (5) 
with a more recent P-L relation fitted by Groenewegen & 
Whitelock(1996, Eq. (2)), 

logP = 2.5/2.59 • log (L/Lq) - 2.71/2.59 . (6) 

Equation (6) is derived for carbon Miras with periods in the 
range 150-520 days. Hence there is still a need to extrapolate 
to longer periods. Winters et al. (1997) point out that these ex- 
trapolated periods are possibly overestimated. The tendency to- 
wards overestimated periods can also be seen in Groenewegen 
et al. (2000, Fig. 3) where (O-rich and) C-rich stars with longer 
periods are measured. These authors point out that both P-L re- 
lations in our discussion are limited to periods of less than 400- 
500 days because they are based on a sample of optically visi- 
ble stars with shorter periods. The extrapolated periods (> 520 



Table 2. Model parameters (L,, T^g and P) of both constant-opacity 
models (prefix "R") and Planck mean models (prefix "P"). Columns 
5 & 6 give the stellar radius, calculated from L, and Ten (Eq. (8)), 
and the gravitational acceleration at the location of the photosphere. 
The last column shows the fraction of the stellar mass contained in 
the model domain of the initial model for constant-opacity models 
and Planck mean models (the latter averaged over the carbon/oxygen 
ratios used), respectively. The stellar mass M, is set to I.OMq in all 
models. The remaining two (out of totally six) model parameters, the 
carbon/oxygen ratio ec/^o and the piston velocity amplitude A Up are 
specified in Tables 4-6 (Sects. 4 & 5.1). 



model 






P 






MJM, 


R-/P- 


[Lq] 


[K] 


[d] 


[%] 




[%] 


07F 


7.0 X 10^ 


2880 


390 


336 


-0.61 


13; 0.16 


lOF 


1.0 X 10" 


2790 


525 


428 


-0.82 


17; 0.16 


lOG 


1.0 X 10" 


2790 


653 


428 


-0.82 


17; 0.16 


13F 


1.3 X 10" 


2700 


650 


521 


-0.99 


20; 0.19 


13G 


1.3 X 10" 


2700 


841 


521 


-0.99 


20; 0.19 



days) given by the P-L relation of Groenewegen & Whitelock 
(1996) should thus be considered an upper limit. 

3. Details of the modeling 

The equations of the RHD3-system are discretized in the 
volume-integrated conservation form on a staggered mesh. The 
spatial discretization of the advection terms is always of sec- 
ond order (van Leer 1977) in this article, contrasting Paper I 
where it was never higher than first order in the drift mod- 
els. The RHD3-system of equations is solved implicitly using a 
Newton-Raphson algorithm where the Jacobian is inverted by 
the Henyey method (cf. Paper 1 for details). 

3.1. Selecting model parameters 

The model is determined by four stellar parameters and two pa- 
rameters defining the piston, there are thus in total six parame- 
ters. The parameters are, the stellar mass M,, the stellar lumi- 
nosity Lt, and the effective temperature T^sf defining the pho- 
tosphere of the initial hydrostatic model. All abundances but 
that of the carbon are assumed solar; the carbon abundance is 
specified through the carbon/oxygen ratio sc/sq. Furthermore 
the period P and the piston velocity amplitude A Up define the 
piston. 

The purpose of this article is to study certain issues con- 
nected to drift that we believe are important for realistic wind 
structures. Therefore we intentionally select combinations of 
model parameters using established parameter relations. The 
sets of model parameters in this study, given in Table 2, are 
mainly selected from the models given in HD97 and HJLA98. 
Note that constant-opacity models are given the prefix "R" in 
the model name, and Planck mean models the prefix "P". In 
selecting model parameters we intentionally leave out the least 
luminous models with the less massive outflows. These models 
are a priori not expected to form winds when adopting drift. 
We restrict the parameter space by not varying the stellar mass 
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M,, which is set to 1 .0 Mq in all models. The piston amplitude 
Awp and the carbon/oxygen ratio sc/sq are currently consid- 
ered "free parameters"; the values we use are specified in the 
model names given in the tables in the results section (Sect. 4). 

In the models presented here it is the dust that initiates and 
drives the stellar wind. The dust formation zone, and hereby 
the wind formation zone, is always inside the model domain. 
An inflow of matter through the inner boundary which is lo- 
cated well below the photosphere (typically at about 0.9/?,) is 
not allowed. One consequence of the lower density present in 
the Planck mean models, compared to constant-opacity mod- 
els, is that a smaller fraction of the stellar mass is located in the 
atmosphere (compare the values in the last column in Table 2). 
Consequently it is currently not possible to model a long-term 
wind evolution since the model domain is depleted quickly of 
material. 

The effective temperature is specified using the radius- 
mass-luminosity relation - derived from an evolutionary model 
-given by Iben(1984)', 

'-'''(w) (Uls) (oool) y 

where the stellar radius R, the stellar luminosity L and the stel- 
lar mass M are given in solar units. The remaining values in 
Eq. (7), the metalicity Z - 0.02 and the ratio of the mixing 
length to the pressure scale height {I /Hp) - 0.90 are taken from 
Bowen & Willson (1991). S is identically unless M > 1.175 
when 5 = 1. The effective temperature is given by, 

CTBT^n- = LJ (AttRI) (8) 

where os is the Stefan-Boltzmann constant. 

The period P is chosen according to the P-L relations men- 
tioned in Sect. 2.3. Model names with an "F" adopt the P-L re- 
lation derived by Feast et al. (1989) (F-models), and those with 
a "G" that derived by Groenewegen & Whitelock (1996) (G- 
models). The relation by Feast et al. is used in all models pre- 
sented in Sect. 4. Results using different periods, represented 
by the G-models, are discussed in Sect. 5.1. 

3.2. Modeling procedure 

The modeling procedure is the same for all models we present 
and is as follows. Each wind model is started from a hydrostatic 
dust-free initial model where the outer boundary is located at 
about 2 Rt. All five dust equations, i.e. the dust equation of mo- 
tion and the four dust moment equations, are switched on at 
the same time. Dust starts to form whereby an outward motion 
of the dust and the gas is initiated. The expansion is followed 
by the grid to about 25 Rt, where the outer boundary is fixed 
allowing outflow. The model evolves for about 50-200 P in the 
case of drift models, and for about 400-900 P in the case of PC 
models. The shorter evolutionary times for the drift models are 
due to the generally larger variations of the dust quantities (in 

' Equation 7 is based on old opacities resulting in overesti- 
mated effective temperatures, the newer relations presented by e.g. 
Wagenhuber & Groenewegen (1998) may be a better choice since they 
are based on more recent opacities. 



Table 3. Model comparison of driftis, models (Cols. 2,3) vs. drift2nd 
models (Cols. 4,5); completing Table 2 in Paper I (where the prop- 
erties for the driftis, models are taken). Only those models that give 
a wind are shown. All wind models show an irregular variability. 
Models that are not run for a longer time are marked with parentheses. 
See Sect. 3.3 for further information. 





drift 


models 


drift^r.,l 


models 


model 


IO'^-<M> 




IO''-<M) 


<"o„) 




[Moyr-'] 


[km s-'] 


[Moyr-'] 


[km s-'] 


A27 


(6.6) 


(40) 


(4.7) 


(40) 


B21 


(11) 


(28) 


(7.5) 


(30) 


B22 


(7.8) 


(28) 


(7.6) 


(32) 


CIS 


15 


21 


12 


23 


CI9 


14 


24 


II 


26 


C20 


13 


25 


II 


28 


DI6 


43 


27 


38 


28 



particular) and the resulting smaller time steps in these models 
(also see Sect. 4.2 in Paper I). One may argue that the terminal 
velocities determined at our outer boundary which is located 
relatively close to the star may be underestimated. However, 
our experience shows that this does not appear to be the case-. 

The massive envelopes of constant-opacity models prevent 
depletion in these models. The density is on the contrary lower, 
and the envelope thereby less massive, in the Planck mean 
models. As a consequence the calculated quantities of the latter 
models are averaged over a time-period shorter than one where 
the depletion becomes relevant (cf. Sect. 4.2). 

3.3. Comparison of numerically improved 

non-pulsating wind models with old models 

In Paper 1 we calculated drift models without stellar pulsations 
using a P' order spatial discretization of the advection terms. In 
this subsection we compare those models (Table 2 in Paper I; 
driftist models) with improved drift models that use a 2"'' order 
advection (drift2nd models). Table 3. Like in Paper I none of the 
improved models simulate stellar pulsations and for compari- 
son the heating term Adrift is not included. Considering the type 
of winds formed in the driftist models, it was found that they 
were all irregular, but so were the PC models that used first or- 
der advection. The new drift2nd models are different in that they 
show much less irregularity and have shorter time-scales in the 
variational patterns. They are, however, all still irregular wind 
models. 

While the measured temporally averaged terminal veloc- 
ities are more or less the same when comparing drift2nd and 
driftist models, the mass loss rates differ The drift2nd model 
A27 has an average mass loss rate that is much lower than that 

- Tests have been performed (with previous versions of the code), 
where the outer boundary of the model is located at larger distances 
from the star. These tests have shown that the uncertainties in the ter- 
minal velocity caused by the location of the outer boundary are small 
compared to the uncertainties in the calculation of the time-averaged 
values. 
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of the driftist model and is about the same as the value of the 
corresponding PC model. All other models have average mass 
loss rates in the range of 72-80%/68-97% of the correspond- 
ing PC/drift models (with a larger uncertainty for the shorter 
evolved B21 and B22 models). 

The conclusion is that with the more accurate discretiza- 
tion of the drift2nd models, compared to the driftist models, 
we find clear changes caused by drift in the wind structure 
of certain models. However, since the properties of both types 
of model winds are of the same order of magnitude and the 
dust still accumulates to the regions behind shocks, there is a 
rough qualitative agreement between driftist and drift2nd mod- 
els. Nevertheless a more accurate discretization is always pre- 
ferred. 

4. Results 

In this section we present and compare the results of new drift- 
and PC models which include the effects of stellar pulsations. A 
discussion of the consequences of the results is given in Sect. 5. 

All quantities given in both tables in this section (and in 
Table 6 in Sect. 5.1) are temporal means calculated at the 
(fixed) outer boundary. To make an estimate of the "degree of 
variability" for each model and quantity we calculate a standard 
deviation, see Sect. 5.2 for a discussion of its relevance and im- 
plications. Note that while all models are calculated for a total 
time-interval which we believe is long enough to find reliable 
average quantities, i.e. 50-400 P depending on the model, we 
do not study long-term variations of the order of several hun- 
dred to thousands of years. Such a study requires total time- 
intervals on the order of thousands of periods and longer 

Depending on the choice of a time-interval used in the 
calculation of the average quantities slightly different values 
emerge. Since a majority of the models in this study are irreg- 
ular it is difficult to determine general conditions for the ex- 
act choice of such a time-interval, instead we choose them as 
long as possible (see, however. Sect. 4.2). We mention at this 
point that the time-interval used for the determination of the 
average drift velocity in our models necessarily is different, cf. 
Sects. 4.1.2 & 4.2.2. 

The two different kinds of gas opacities we apply give in- 
trinsically different wind structures and we therefore present 
the results of each of these two models in separate subsections. 
We first study the results of constant-opacity models and there- 
after the results of (and a comparison with) the Planck mean 
models. To illustrate differences we have chosen to put results 
of different nature in the same figure. To aid the further dis- 
cussion we here give the following symbol key for these fig- 
ures (Figs. 1-3, 6, 7): PC models are always represented by 
open symbols, and drift models by filled symbols. Moreover 
constant-opacity models are represented by circular symbols, 
stars, and downwards pointing triangles (o, ©, @, •, r), and 
Planck mean models by rectangular symbols and upwards 
pointing triangles (□, a, m, a). Planck mean drift models are 
plotted in gray to distinguish them more clearly from constant- 
opacity drift models. Triangular symbols represent models with 
a varied period P. Further details are given in the respective ta- 
bles and figures. 



4.1. Results of constant-opacity models 

The results of the constant-opacity models are given in Table 4. 
PC models are given in the upper part and drift models in the 
lower part; note that we adopt the same model name for PC 
models and drift models with the same model parameters. The 
standard deviation os of the average values is specified for the 
first five quantities (but not for the drift velocity, cf. Sect. 4.1.2). 
We give the actual dust/gas density ratio in Col. 8 (instead of 
calculating it from the degree of condensation like in HD97 & 
HJLA98, which only works correctly for PC models). 

4.1 .1 . Comparison of new PC model results with 
previous results 

Before we study the drift models we comment on the agree- 
ment of the values given for the PC models in Table 4 (with an 
'H' in the last column) with the corresponding values given in 
HD97. The physics of the part of the RHD3-system required to 
run PC models is not modified, and the result should therefore 
be the same (or very similar). The PC model values shown in 
bold face in Table 4 are values where the difference is > 10%. 
Comparing old and new values we find that the agreement gen- 
erally is fine. 

A factor that explains several of the differing values is 
that all new PC models are calculated for significantly longer 
time-intervals than the models in HD97. Models that show dif- 
ferent average values due to the longer time-interval include 
R07FU2C20, R10FU2C15, R10FU2C16 and R13FU2C14. 
The latter two models show large variations during parts of the 
evolution and the last model relaxes to a stable pattern of vari- 
ability only after several hundred piston periods. One model 
where the above arguments do not apply and the average val- 
ues of the new models are different is R10FU2C20. For this 
model all average values are about ten percent larger, except 
for the mass loss rate which is 27% larger In this case it is the 
selection of the time-interval used in the calculation of the av- 
erage properties that is important (compare with the situation 
of the Planck mean models in Sect. 4.2.1). 

4.1 .2. Comparison of drift models with PC models 

A direct comparison of the values of PC models and the cor- 
responding values of drift models in Table 4 shows significant, 
but not dramatic changes. In this subsection we study correla- 
tions we can see by comparing values of drift and PC models. 

The direct comparison between PC models and drift models 
is aided by a plot of the ratios of the values. Fig. 1. In this 
subsection we only discuss the models represented by circular 
symbols and stars (o, ©, @, •, ★) in Figs. 1-3. 

Comparing the model ratios in the panels of Fig. 1 we 
find ratios both significantly smaller and larger than 1. With 
an exception of the ratios of two models that seem to be sep- 
arated from the rest of the model ratios in Figs, la-c (see be- 
low), the other models are grouped rather tightly. The ratio of 
the mass loss rates is slightly lower than 1 for the group of 
constant-opacity models, and the terminal velocity ratio is at a 
constant level close to 1 . The degree of condensation is fairly 
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Table 4. Model quantities averaged at the outer boundary for constant-opacity models (cf. Sect. 4.1). The models are named by adding a suffix 
to the respective model name in Table 2. The suffix is a combination of two model parameters, the piston velocity amplitude A Up (Un, n in 
kms"') and the carbon/oxygen ratio ec/sq {Cnn, nn is 10 ■ bc/eq)- The last but two column, r,o,, gives the total time-interval of each model 
calculation. The numbers associated with models run for a shorter time-interval are less reliable than the others as the means are taken over 
short time intervals. The last but one column gives the type of wind: i, irregular wind; Ip, periodic wind; /q, quasi-periodic wind; t, transition 
model; — , no wind. / (e N) shows the (multi-)periodicity of dust shell formation in the unit of the piston period P. A tilde (e.g. '2q') in the 
last but one column indicates a correspondence with the characterization only during a part of the calculated time-interval. The values shown 
in bold face for the PC models indicate that the value differs from the corresponding value given by HD97 (by > 10%). The symbols in the 
last column show whether the respective drift models show an increased/decreased (a/v) mass loss rate, or new wind (®) when compared to 
the corresponding PC model. For the PC models an 'H' indicates that the same model parameters were used in a model in HD97. The symbol 
printed in subscript in the last column indicates how the respective (differing) model is illustrated in Figs. 2, 3, 6 & 7. 
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Fig. 1. Ratios of four of the drift/PC model quantities given in Tables 4 & 5 plotted as a function of the respective PC model quantity. From 
the left the averaged quantities are: a) the mass loss rate {M)\ b) the terminal velocity {u^); c) the degree of condensation (/cond); and d) the 
dust/gas density ratio (pilp). Note that all plots are logarithmic on the y-axis and of the same scale. Only models where both PC models and 
drift models have formed a wind are shown. The ratios of constant-opacity models are represented by open circles o, and Planck mean models 
by filled squares ■, respectively. The two constant-opacity models showing an increased mass loss rate when allowing drift are indicated in 
all panels with the symbol @. The region enclosed by the horizontal lines, 0.9 < ratio < 1.1, indicates an insignificant decrease/increase of the 
respective ratio. Note that each set of constant-opacity models and Planck mean models tend to form groups; excepting the two models. 

Fig. 2. Model average quantities as a function of the mass loss rate (M). The panels show: a) & b) the degree of condensation (/cond) and 
c) & d) the optical depth of the dust (tj). The panels on the left show PC models (open symbols) and the panels on the right drift models 
(filled symbols). Furthermore, constant-opacity models are represented by circles, stars, and downwards pointing triangles (o, ©, @, •, t), 
and Planck mean models by squares and upwards pointing triangles (□, ffl, ■, k). A key to the symbols is given just before Sect. 4.1, and in 
Tables 4-6. Note that the PC models without a corresponding drift model (indicated by © & ffl) all have low mass loss rates. The drift models 
are as a group better correlated with the mass loss rate in the degree of condensation (b) than the PC models are (a). 
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constant at about fcond - 0.7 for PC models; the drift mod- 
els show a larger scatter of values, where the individual values 
always are smaller than those of PC models. The two mod- 
els that are separated from the group, i.e. models R10FU2C15 
and R13FU2C14, both share an increased drift/PC model ra- 
tio in the mass loss rate, the terminal velocity, and the dust/gas 
density ratio. In addition both models share a lower degree of 
condensation and a low terminal velocity compared to the other 
constant-opacity models. These two models are discussed fur- 
ther in Sect. 5.2.1. 

To illustrate correlations between the average model quan- 
tities we show Figs. 2 & 3. Figure 2 relates the degree of con- 
densation (/cond) (upper panels) and the dust optical depth (rd) 
(lower panels) to the mass loss rate. The degree of condensa- 
tion measures how easily dust forms and the dust opacity shows 
what influence the amount of dust in the wind has on the wind 
structure. In the current context we mention the three models 
that give a higher mass loss rate in the drift models (i.e. mod- 
els R10FU2C15, R13FU2C13 & R13FU2C14), indicated with 
stars (★, drift models) and open rings (@, PC models). 

Figure 2 shows that for the PC models there is no corre- 
lation of the degree of condensation with the mass loss rate, in 
addition three models (R10FU2C15, R13FU2C14 and the tran- 
sition model R07FU2C18) share a significantly lower degree 
of condensation. The drift models, however, seem to follow a 
trend of increased mass loss rate with the degree of condensa- 
tion (or vice versa). Both the PC models and the drift models 
show a correlation of the mass loss rate with the dust optical 
depth. 

As pointed out by HD97 there should be a close correla- 
tion in dust-driven winds between the terminal velocity and the 
strength of the radiation pressure relative to the gravitational 
pull. In Fig. 3 we plot the terminal velocity (lower panels) and 
the mass loss rate (upper panels) against the quantity (a). We 
adopt the same expression for (a) as used by IID97 (Eq. (1)), 
where it is defined at the outer boundary; (a) is proportional 
to the ratio of the dust radiation pressure term (/lad.d) and the 
gravitational force acting on the gas (/grav.g) in the (PC model) 
equation of motion (see Eq. (1) in Paper I). The terminal ve- 
locity is, however, affected in the entire wind formation region 
and a more certain result might be given with a radially aver- 
aged (a). WOO calculate a similar property (at) for the inner- 
most dust shell using the radiative pressure on both the gas and 
the dust. Note that (a), by definition, is not a dimensionless 
quantity, and therefore a value > 1 (< 1) does not indicate that 
the radiation pressure (gravitational pull) dominates. 

Figure 3 shows that the mass loss rate is not correlated to 
{&}, neither for the drift models nor the PC models, also com- 
pare the scatter of models in Fig. 3a with Fig. 2 in WOO. A 
strong correlation of the velocity to (a) is seen in the PC model 
plot of the terminal velocity in Fig. 3c, in agreement with Fig. 3 
in HD97. The agreement with Fig. 3 in WOO is less obvious. 
Because of the larger scatter of values the drift models are less 
well correlated in both panels; which could be an indicator of 
the weaker coupling between gas and dust in these models. 
Note that the three models that give an increased mass loss rate 
in drift models all are situated in the region of low {a}, for both 
PC and drift models. 



The drift velocity in Col. 12 in Table 4 requires a separate 
discussion. Since the dust velocity often shows large and dra- 
matic variations in regions in front of shocks (see discussion 
in Sect. 4.2 in Paper I), it is not possible to calculate an av- 
erage value for all models using the same method as for the 
other quantities. Instead we base this average on the longest 
time-interval that does not show the largest unphysical values. 
A small average drift velocity indicates a frequent passing of 
shocks (where the drift velocity is low, cf. Sect. 4.2 in Paper I), 
and a large value correspondingly few shocks. The values in 
the table show a tendency towards larger drift velocities with 
less massive winds. We do not give a standard deviation for 
this quantity. 

4.2. Results of Planck mean models 

The results of the Planck mean models are given in Table 5. 
The PC models are presented in the upper part and drift mod- 
els in the lower part. In contrast to constant-opacity models, 
Planck mean models have a lot less mass in the model domain 
and are more rapidly depleted of material. All PC models in 
this subsection are evolved for a period long enough to show a 
significant depletion. Consequently it is not suitable to use the 
same long time-intervals in the calculations of averages that 
we used for constant-opacity models, instead we set the upper 
limit of the time-interval at a time before the depletion becomes 
significant. 

4.2.1 . Comparison of new PC model results with 
previous results 

Like for the constant-opacity models in Sect. 4.1.1 we com- 
ment on the agreement of the values given for the PC models 
(indicated with an 'H' in the last column) in Table 5 with the 
corresponding values given in HJLA98. The result should be 
the same (or very similar). The PC model values shown in bold 
face in Table 5 are values where the difference is > 10%. 

The length of the time-interval used in the calculation of 
average quantities is a factor that can explain several of the 
different values in Sect. 4.1.1. Because of the fast depletion of 
the model domain in the Planck mean models this it is not the 
cause of discrepancy here. Instead we point out that it is the 
detailed selection of the time-interval used in the calculation of 
the averages that is crucial. We select this upper time-interval 
limit for models showing depletion to be where the average 
decrease in the mass loss rate is insignificant. This point can be 
very difficult to define in winds of irregular variability and it 
is therefore wise to choose a shorter period. The approximate 
time (in piston periods) printed in subscript in the last but two 
columns in Table 5 indicates this upper limit. 

All models that differ by > 10% in Table 5, except 
P07FU6C14, are found to differ due to the selection of the time 
interval. Model P07FU6C14 evolves very slowly to a stable 
variability pattern, and it could be that we measure the quan- 
tities during different states of the wind (using different time- 
intervals). 



C. Sandin & S. Hofner: Three-component modeling of AGB star winds. II 



9 



Fig. 3. Model average quantities as a function of (a), a quantity characterizing the strength of the radiation pressure relative to the gravitation. 
The panels show: a) & b) the mass loss rate (M> and c) & d) the terminal velocity (Uoo). The panels on the left show PC models (open symbols) 
and the panels on the right drift models (filled symbols). Furthermore, constant-opacity models are represented by circles, stars, and downwards 
pointing triangles (o, ®, ©,•,★, r), and Planck mean models by squares and upwards pointing triangles (□, ffi, ■, a). A key to the symbols is 
given just before Sect. 4. 1 , and in Tables 4-6. Note that the PC models without a corresponding drift model (indicated by © & ffl) all are grouped 
in the lower-left corners in a & c. The Planck mean and constant-opacity PC models are clearly more separated as groups compared to the 
corresponding (opacity) groups of the drift models. Also compare c with Fig. 3 in HD97. 

Table 5. Model quantities temporally averaged at the outer boundary for Planck mean models. The numbers printed in subscript in the T,at 
column indicates the upper limit in the time-interval used in the calculation of the average quantities of the corresponding model (cf. Sect. 4.2. 1). 
The values shown in boldface for the PC models indicate that the value differ from the corresponding value given by HJLA98 (by > 10%). For 
the PC models an 'H' in the last column indicates that the same model parameters were used in a model in HILA98. See the caption of Table 4 
for further details. 
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4.2.2. Comparison of drift models witli PC models 

Due to the overall lower density in the Planck mean models 
these models form a lot less dust than the constant-opacity 
models do. The results of this are for example a lower aver- 
age degree of condensation and a lower average mass loss rate 
in the former type of models (cf. HJLA98 and Helling et al. 
2000). A first difference seen in the drift models of Table 5 is 
that only three of the sets of HJLA98-PC model parameters are 
found to form a drift model wind. None of the PC models with 
a mass loss rate < 1 .0 x lO-*" M0 yr-' successfully forms a drift 
model wind. With fewer drift models the statistics is on the one 
hand less accurate and correlations become less certain. On the 
other hand, with the current selection of model parameters, we 
get stronger conditions on the limit for wind formation in drift 
models (cf. Sect. 5.2). 



The results of all Planck mean models are plotted together 
with the constant-opacity models in Figs. 1-3. In this subsec- 
tion we study models represented by the rectangular symbols 
(□, ffl, ■) and compare with the relations of the constant-opacity 
models. We again find both increased and decreased values in 
Fig. 1 . There is, however, no model (or group of models) show- 
ing a significantly increased mass loss rate (like two constant- 
opacity models do). While the values of the PC models group 
clearly around 13 kms-' in Fig. lb, around 0.2 in Fig. Ic, and 
around 0.8x10-^ in Fig. Id, there is a larger scatter in the values 
of the drift models. There is no evident correspondence with the 
scatter of the constant-opacity models. Like for the constant- 
opacity models the (gray) triangles (of the Planck mean models 
in Table 6, Sect. 5.1) show a distribution similar to the rest of 
the drift models. 
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In Fig. 2 (and Fig. 3) it seems that the Planck mean PC mod- 
els are scattered over a larger region than the constant-opacity 
PC models are. However, if those models that do not give a 
corresponding drift model are excluded (indicated by a), six 
tightly grouped models are left. For the drift models it appears 
that the Planck mean models correlate well with the constant- 
opacity models, in particular in the dust optical depth (Fig. 2d). 
The drift models in Fig. 3 share the properties of - and show 
a distribution similar to that of - the constant-opacity models. 
Excepting one model, P10FU6C18, the Planck mean models 
compared to the constant-opacity models, on average share a 
lower value of (a) (this is caused by a lower dust/gas density 
ratio. Fig. Id). 

Note that all PC models with a terminal velocity below 
lOkms ' either lack a corresponding drift model, or show sig- 
nificantly different properties when drift is included. The av- 
erage drift velocities in the six drift models are all fairly low, 
and two models are even multi-periodic during parts of the cal- 
culated time-interval (cf. the discussion in Sect. 5.1). In addi- 
tion to the variability associated with the shell structure of the 
winds, the dust velocity shows small-amplitude variations at 
the outer boundary for all drift models in Table 5 (see Sect. 5.2). 

5. Discussion 

In the current models drift does not change the fundamental 
processes behind the wind formation mechanism, but is still 
able to alter the conditions for dust formation. In wind models 
allowing drift the dust tends to accumulate to dense regions be- 
hind shocks. The reason for this relocation is that the gas-dust 
collisional interaction is weak in inter-shock regions and the 
dust, accelerated by the radiation pressure, drifts with respect 
to the gas until the gas-dust interaction is strong enough to slow 
the drift down. The dense regions behind shocks provide such 
conditions (cf. Paper I). 

Figure 4 shows the spatial structure of the constant-opacity 
drift and PC models R10FU2C18. The quantities are selected to 
illustrate the difference in the spatial wind structures of the two 
models. The accumulation of dust to the regions behind shocks 
is evident in the degree of condensation and the dust/gas num- 
ber density ratio; the dust density is several orders of magnitude 
lower in the inter-shock region, where the drift velocity on the 
contrary is higher. The shock/inter-shock gas density ratio is 
much lower than that of the dust. Although the dust densities 
of the two models diff'er, the gas structures appear similar in the 
gas velocity and the gas density. 

Figures 4b,f show that both the average grain radius and the 
grain abundance vary more strongly in the drift model. The av- 
erage grain radius is typically below 10"^' cm, justifying the use 
of the small particle limit in the calculation of the dust opaci- 
ties. The high values in the innermost peak in the plot of the 
average grain radius, at 2 are irrelevant since the degree of 
condensation, and therefore the dust opacity, at this location is 
negligible (compare Fig. 4b with Fig. 4d). 

In most cases of our study the consequences of the accu- 
mulation of dust to shocks tends to be (on average) narrower 
regions of efficient nucleation and grain growth. This in turn 
leads to smaller total amounts of formed dust, a less efficient 



wind formation, and lower mass loss rates compared to winds 
formed in PC models. However, this is not necessarily the case 
if the dust formation in a certain PC model is inefficient. In a 
few cases the accumulation of dust to narrow regions (or shells) 
in drift models seems to provide more favorable conditions for 
the wind formation with a resulting increased mass loss rate. 

In this section we study four topics of the wind formation 
where the influence of drift is of a primary concern. We first 
study the influence of the piston period in the next subsection. 
In Sect. 5.2 we measure and compare the variability of the wind 
models. The comparison with a derived mass loss rate fit for- 
mula is discussed in Sect. 5.3, both for the PC and the drift 
models and for constant-opacity models vs. Planck mean mod- 
els. Finally we comment on the importance of the heating by 
drift for the energy balance in drift models in Sect. 5.4. 

5.1. Influence of the piston period 

We believe it is important to follow up the study regarding 
the influence of the piston period on the wind structure car- 
ried out by WOO (Sect. 3.5.6) and HD97 for the case of drift 
models. In particular regarding Planck mean models, for which 
there is no similar study. All models discussed this far are 
calculated with a P-L relation derived for M-type Miras by 
Feast et al. (1989, we refer to these models as F-models). 
A more recent relation derived for C-rich Miras, derived by 
Groenewegen & Whitelock (1996, G-models), reveal longer 
periods (see Sect. 2.3). The results of seven G-models are pre- 
sented in Table 6 together with the results of the respective 
F-models for comparison. The wind properties that differ by 
more than 10% from the corresponding F-model are indicated 
in boldface. The piston periods of the (F- and) G-models are 
given in Table 2. 

In the constant-opacity models (upper half in Table 6) we 
note that the values of the G-models primarily differ in the mass 
loss rate and the degree of condensation. Quantitatively the 
changes are of the same order of magnitude as found by HD97. 
The exception is model R10GU2C18 where the degree of con- 
densation is about half of that of R10FU2C18. Two G-models, 
R10GU2C18 and R10GU2C15 also show multi-periodic vari- 
ations where the F-models do not. These two models and the 
F-model R13FU2C13 are the only constant-opacity drift mod- 
els in this study that are multi-periodic. In the latter model a 
dust shell forms every third piston period. 

Planck mean models (lower half in Table 6) are more sen- 
sitive to the period than the constant-opacity models are. The 
values of all three G-models differ significantly compared to 
the corresponding F-models. One model, P13GU6C16 is even 
multi -periodic during parts of the evolution. This model and the 
two F-models P10FU6C18 & P13FU4C16 are the only three 
Planck mean models (among both drift and PC models) where 
we find multi-periodicity. The changes show both increasing 
and decreasing values in all quantities. Additionally in the mass 
loss rate, the terminal velocity, and the dust optical depth all G- 
models show higher values. 

We conclude that while constant-opacity models are only 
httle affected by the piston period in some models (in accor- 
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Fig. 4. Radial structure of the inner parts of the wind for the drift (solid line) and PC (dash-dot-dotted line) models R10FU2C18. The panels 
show: a) the gas velocity u; b) the average grain radius (r^); c) the gas density p; d) the degree of condensation f^oud', e) the drift velocity vd; and 
f) the grain/gas number density ratio na/n,. The accumulation of the dust to regions behind shocks in drift models is seen in e.g. the dust/gas 
number density ratio (f) and the degree of condensation (d). Note the larger variations in the dust quantities of the drift model (right panels). 

Table 6. Results of drift models using two different P-L relations. F-Models are indicated with an "F" in the model name (first column), and 
G-models with a "G", cf. Sects. 3.1 & 5.1. The symbols in the last column indicate whether the mass loss rates of the G-models have an 
increased/decreased (A,v) value compared to the corresponding F-models. The values shown in boldface for the G-models indicate that the 
value differ from the corresponding F-model value (by > 10%). The symbol printed in subscript for the F-models in the last column indicates 
how the respective model is illustrated in Figs. 2, 3 & 6. See the caption of Table 4 for further details. 
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dance with the results of WOO and HD97), Planck mean mod- 
els are sensitive to the period to a larger extent with significant 
changes in all wind properties. 

5.2. The importance of variability in the wind 

The temporal evolution at the outer boundary is shown for 
the drift and PC models R13FU4C14 (constant-opacity mod- 
els, left panels) and P13FU6C14 (Planck mean models, right 
panels) in Fig. 5. All corresponding panels on the left and the 
right are drawn at the same scale on the axes for easy compar- 
ison. The same model parameters are used in both the Planck 
mean models and the constant-opacity models, with the excep- 
tion of the piston amplitude which is higher in the two Planck 
mean models. The higher value is used as a partial compensa- 
tion for the lower densities in Planck mean models. The result 
of the less efficient dust formation in the Planck mean model 
P13FU6C14 is clearly evident through all panels in the figure. 

The large variability of the drift models make a direct com- 
parison with the PC models difficult in terms of which wind has 
larger average values. Although the variations in the terminal 
velocity and the mass loss rate seem to be of the same magni- 
tude and time scale for the models R13FU4C14, they differ in 
both the PC and drift models P13FU6C14. Note that the varia- 



tional patterns in the dust quantities of the drift models (lower- 
most four panels for both pairs of models in Fig. 5) differ signif- 
icantly from the corresponding patterns of the PC models. For 
the two drift models shown in this figure less dust is formed, 
the grains are smaller, and the mass loss rate is decreased in 
comparison with the values of the respective PC model. 

In the lowermost panels in Fig. 5 we see periodic small- 
amplitude variations superposed on the drift velocity. These 
variations do not appear in the gas velocity and last for the en- 
tire calculated time-interval. They are larger in the less dense 
inter-shock regions. In addition, compared to the constant- 
opacity model, they are more pronounced in the less massive 
wind of the Planck mean model. The period of these varia- 
tions is for both models equal to the piston period. By the 
strong coupling to the radiation field the dust senses the vari- 
ations in the luminosity caused by the radial oscillations of 
the piston. The small-amplitude variations remain in the dust 
velocity as the coUisional coupling between the gas and the 
dust is not strong enough to flatten them out. These variations 
ought to be present in all models, but since the long-term vari- 
ations in the drift velocity are larger in several models they 
are not always obvious. Corresponding modulations are visi- 
ble in all Planck mean models except the multi-periodic model 
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Fig. 5. The temporal evolution at the outer boundary. The panels on the l.h.s. show the constant-opacity drift and PC models R13FU4C14 
(solid and dash-dot-dotted lines, respectively), and the panels on the rh.s. the Planck mean drift and PC models P13FU6C14 (solid and dash- 
dot-dotted lines, respectively). With the exception of the piston amplitude Awp, which is larger in the Planclc mean models, all other model 
parameters are the same. The panels show (from the top): the terminal velocity the mass loss rate M; the degree of condensation f^aM, the 
dust/gas density ratio p^/p; the mean grain radius (rj); and the drift velocity v^. Note the much larger variations present in both drift models 
compared to the PC models, also compare with Fig. 6. 

Fig. 6. Temporal variations of average properties {q) measured by the fluctuation amplitude r = ajq (Sect. 5.2). From the left the figure shows: 
a) the mass loss rate (M); b) the terminal velocity c) the degree of condensation (/cond); and d) the dust optical depth (tj). All plots are 
drawn to the same scale. Filled symbols represent drift models (lower panels) and open symbols PC models (upper panels), constant-opacity 
models are represented by circles, stars, and downwards pointing triangles (o, ©,©,•,*, ▼), and Planck mean models by squares and upwards 
pointing triangles (□, ffl, ■, A). Furthermore, PC models without a corresponding drift model for the same set of model parameters are indicated 
with the symbol ® or ffl. Models for which the drift models have a larger average mass loss rate are indicated with open rings (@; PC models), 
and stars (★; drift models). Triangles represent G-models (cf. Sect. 5.1). Note the distinctly smaller scattering and larger variations of the values 
of the drift models compared to the values of the PC models. 



P10FU6C18. In addition they appear in three constant-opacity 
models: R10GU2C15, R13FU4C14 and R13FU2C16. 

To quantify the variability of the model properties we have 
calculated an (absolute) standard deviation (cr^) for each aver- 
age outflow quantity (except the drift velocity, see Sect. 4.1.2) 
given in Tables 4-6. However, a relative fluctuation amplitude 
is preferred in comparisons between different models, see be- 
low. The relative fluctuation amplitude is simply the standard 
deviation divided by the average value (r - cr^/q, i.e. a "rel- 
ative error"). A large variability, typical of a time-dependent 
outflow, is represented by a large fluctuation amplitude. And 
a small variability, typical of a "stationary" outflow, by a small 
fluctuation amplitude. In a similar argument WOO used the stan- 
dard deviation to quantify the stability of the wind. 

In the time sequences shown in Fig. 5 we saw a tendency 
of a larger variability in the two drift models, compared to the 
corresponding PC models. The fluctuation amplitude is shown 
in Fig. 6 as a function of four of the average properties of 
Tables 4-6. With Fig. 6 we see that the variations generally are 
slightly larger to significantly larger in the drift models of most 
drift/PC model pairs. 

By marking the PC models that lack a corresponding drift 
model (with the symbols ffi & ffl) we find that these models, 
compared to the other models, either show smaller variations 
or lower average values. Apparently they share several proper- 
ties. For instance, all PC models in this group share a terminal 
velocity < 10km s ', no other PC model does. Furthermore, 
these models tend to have a mass loss rate in the lower end of 
the range. We point out that we do not attempt to model winds 
where we initially are certain that there is no drift model wind. 
The statistics of the subset of PC models without corresponding 
drift models is therefore not sufficient for generalized quantita- 
tive conclusions in this respect. 

We first study the constant-opacity models in the following 
subsection and then Planck mean models in Sect. 5.2.2. 

5.2.1 . Variability of constant-opacity models 

In this subsection we discuss constant-opacity models shown 
with circular symbols, stars, and downwards pointing black tri- 



angles (o, ©, @, t) in Fig. 6. The basic properties of the 
models are studied in Sect. 4.1.2. 

The variability of the drift models are for most cases about 
the same or larger than the variability in the PC models in 
Fig. 6. Two PC models show variations in the terminal veloc- 
ity that are larger than those of the bulk of the drift models. 
These models, R07FU2C18 (0) and R13FU2C14 (@), together 
with R10FU2C15 (@) all have a lower terminal velocity and 
show a lower variability in tj compared to the rest of the PC 
models. Moreover, the latter two models, showing an increased 
mass loss rate in the corresponding drift models, both have a 
degree of condensation below /cond - 0.5. Model R07FU2C18 
barely forms a wind using drift (indicated by a 't' in Table 4). In 
contrast the two corresponding drift models R10FU2C15 and 
R13FU2C14 correlate well in the variability with the other drift 
models. 

To understand exactly what model parameters result in drift 
models with higher mass loss rates we would need a large sam- 
ple of models that do not show any corresponding drift model. 
It would with this sample of models be possible to separate the 
PC models without drift model winds from the corresponding 
models giving a higher mass loss rate in drift models. These 
latter models can in turn, using the arguments in the previous 
paragraph, be separated from the rest of the models that give a 
similar or lower mass loss rate in the drift models. 

In this context we mention the drift model R10FU2C13, 
which lacks a corresponding PC model wind. One attempt to 
form winds with a lower carbon/oxygen ratio, sq/so - 1.25, 
in a drift model with the otherwise same parameters as model 
R13FU2C13 failed. Similarly a second attempt with ec/eo = 
1 .45 in a (PC & drift) model with the otherwise same param- 
eters as model R10FU2C15 also failed. The properties of the 
model with the closest model parameters, with a higher car- 
bon/oxygen ratio, R13FU2C14 shows that the added freedom 
in the drift model is crucial in forming a wind with these pa- 
rameters. 

5.2.2. Variability of Planck mean models 

In this subsection we discuss Planck mean models, shown with 
rectangular symbols and upwards pointing triangles (□, ffl, ■, a) 
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Fig. 7. Ratios of the average mass loss rates (M) to the mass loss fit 
Mfit given in Eq. (9) as a function of: a) & c) the effective tempera- 
ture Jeff; and b) & d) the stellar luminosity L. The upper panels show 
PC model ratios and the lower panels additionally drift model ratios 
(for comparison). The plot symbols are the same as in Fig. 6. The 
horizontal lines indicate the mean ratio for each type of model wind. 
This figure (specifically the shaded region) should be compared with 
Figs. lb,c in Wachter et al. (2002). The Planck mean models are drawn 
in light gray symbols in the upper panels to emphasize that Eq. (9) is 
derived for constant-opacity PC models; cf. Sect. 5.3. 

in Fig. 6. The basic properties of these models are studied in 
Sect. 4.2.2. A comparison of the distribution of model values 
in the figure between Planck mean models and constant-opacity 
models show that the scatter is larger for the former type of 
models, both for PC and drift models. 

All Planck mean PC models that do not have a correspond- 
ing drift model (ffl) show a terminal velocity (mm) S lOkms"'. 
In addition the same four models show low mass loss rates and 
low variations in all quantities. The variations are particularly 
low in the degree of condensation where no other but these four 
Planck mean models show a fluctuation amplitude r < 0.05. 
This indicates the importance of a time-dependent wind struc- 
ture for forming drift models. We conclude that Planck mean 
PC models of a low terminal velocity and low variations in the 
degree of condensation are the models that are least likely to 
form a wind in a drift model. 

5.3. Comparison of the mass loss rates with a mass 
loss rate fit formula 

The calculation of wind models is a time consuming process. 
And it may in some cases be preferable if the calculation of 
the mass loss rate can be replaced by a function with a small 
number of parameters. Wachter et al. (2002) derived a mass 
loss rate fit formula based on the wind models in WOO, that are 
similar to our constant-opacity PC models, and the P-L relation 
of Groenewegen & Whitelock (1996) (see Sect. 2.3). Their ex- 
pression allows the calculation of the mass loss rate using only 
four stellar parameters (the pulsation period P is below indi- 
rectly included in the L-term through the P-L relation), 

logMfit = -4.52- 6.81 log(reff/2600K) 

-1-2.47 log(L/10'*Lo) - 1.951og(M/Mo). (9) 

The units of Mfjt are Mq yr Our models are well within the 
range of stellar parameters covered by Eq. (9). In this subsec- 
tion we discuss how the output of our models compares to this 
expression. 

The mass loss rate ratio (M)/Mfit is illustrated for the mod- 
els presented in this article in Fig. 7. A comparison of the 
constant-opacity PC models in this figure (illustrated by open 
circular symbols (o, ©, @) in the upper panels) with Figs. lb,c 
in Wachter et al. (2002) reveals systematically lower mass loss 
rates in our models. This discrepancy is most likely caused by a 
smaller levitation in our models. The models used by Wachter 
et al. have a piston velocity amplitude Awp = 5kms"' com- 
pared to 2 km s 'used in most of our constant-opacity models. 



Fig. 8. Radial structure of the Planck mean model P10FU6C18 (solid 
line) and the constant-opacity model R10FU2C18 (dash-dot-dotted 
line). The shown quantities are: a) the gas velocity; b) the gas den- 
sity; c) the drift velocity; d) the drag force; e) & f) the energy ratios 
(see Sect. 5.4). This figure illustrates the importance of the heating by 
drift (gdiift) for the energy balance of the respective model. The heating 
by drift is negligible in most parts of the Planck mean model (lower- 
most two panels). On the contrary it is comparable in magnitude to the 
other energy source terms in the constant-opacity model. Even so we 
do not find any significant effect on the wind structures caused by this 
term, cf. Sect. 5.4. 

The two pairs of models with 2 & 4kms ' but otherwise iden- 
tical parameters (see Table 4) demonstrate the influence of A 
- and consequently the levitation - on the mass loss rate. 

The densities in the Planck mean models are much 
lower compared to the corresponding values in the constant- 
opacity models. And the resulting mass loss rates are much 
lower (Sect. 4.2 & HJLA98; illustrated by squares (n,m) in 
Fig. 7). The mass loss rate also differs when drift is allowed 
(Sects. 4.1.2, 4.2.2; lower panels in Fig. 7), but the change is 
less obvious in the figure above. The mass loss rate fit in its 
present form is inapplicable to these winds as it does not ac- 
count for these two effects. However, the computational effort 
needed to calculate a new relation is hardly justified since fu- 
ture further physically improved models possibly quickly will 
invalidate such a relation. 

5.4. The influence of heating by drift 

The importance of heating by drift (^drift) for the wind struc- 
ture was studied by KGS94, who found the heating term (cf. 
Sect. 2. 1) to be an important factor in the wind energy balance 
in the context of the stationary winds they modeled. They re- 
ferred to the term as qfnc- In the models of the first article in 
this series. Paper I, we did not include this term as its inclusion 
is connected to several numerical difficulties. For comparison 
with those models it is not included in the models discussed in 
Sect. 3.3 in this article, but it is included in all other models. 

Places in the wind where the heating by drift could be ex- 
pected to be important are in the dense regions behind shocks 
where the drag force is strong. However, the drift velocity is 
simultaneously low at these locations, and it is therefore diffi- 
cult to make a simple estimate of the importance of this heat- 
ing term compared to the other terms in the energy equation. 
The nonlinear character of the terms in the system of equations 
describing the wind in connection with a highly variable struc- 
ture, makes it difficult to predict when and where in the wind 
Adrift is important. 

With Fig. 8 we illustrate the complex influence of the 
heating by drift. The radial plot shows the Planck mean 
model P10FU6C18 (soUd line) and the constant-opacity model 
R10FU2C18 (dash-dot-dotted line). The uppermost two pan- 
els in Fig. 8 depicting the gas velocity u and the gas density p 
show the locations of the shocks. The middle panels depicting 
the drift velocity vd and the drag force /d,ae are plotted with 
the same scale and illustrate the importance of the drag force 
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in relation to the drift velocity in the Adrift term. The lowermost 
two panels show the relative strength of the heating by drift 
Adrift compared to the term describing the work done on the gas 
by the velocity field (PV-m, Fig. 8e), and the energy exchange 
with the radiation field {AnpK„{J - Sg), Fig. 8f; where is the 
gas source function) respectively. In both these panels a ratio 
less than one indicates that the heating by drift is smaller than 
the respective heating term. Note the logarithmic scale in both 
these panels. 

The spatial structures in each model vary significantly with 
time and the presented instants do not represent all possible 
states of the plotted quantities. They do, however, give an im- 
pression of the typical features and range of values. Figs. 8c,d 
show that the variation of the drag force is more important to 
the q^yift term than the drift velocity is. The lowermost two pan- 
els show that the same term has a negligible influence on the 
energy balance in the Planck mean model; the energy ratios are 
only close to one in the shocked regions. The situation is ob- 
viously different in the denser constant-opacity model, where 
the heating by drift seems to be relevant for larger zones of the 
model. Nevertheless, regarding the average wind properties we 
do not find any measurable "significant" differences between 
models calculated with and without Adrift, neither in any of the 
Planck mean models nor in any of the constant-opacity models 
presented in the tables in Sect. 4. 

We conclude that although the heating by drift seems to be 
of a minor importance for the models we have studied here, it 
can (possibly) not safely be neglected in general. The term Adrift 
should therefore be included in the calculations. 

6. Conclusions 

In this paper we have studied the effects of drift between dust 
and gas in time-dependent wind models of cool long-period 
variables. The drift and its consequences for the wind struc- 
ture have been discussed for several decades. However, until re- 
cently these discussions were almost exclusively based on sta- 
tionary wind models. These models, and the associated struc- 
tures, showing a smooth distribution of gas and dust may be 
misleading when trying to estimate the effects of drift, as stud- 
ies using time-dependent models have demonstrated within the 
last few years. In addition, among the existing dynamical wind 
models only a few actually describe the formation of the stellar 
wind using a self-consistent description of the gas, the dust and 
the radiation field. In this article we have presented models us- 
ing such a description (i.e. the approach introduced by Sandin 
& Hofner 2003, Paper I). 

With this study we wanted to find out how the wind proper- 
ties, and the variability of the wind structure, change as a con- 
sequence of including drift in the calculations. Compared to 
Paper I we now include the effects of stellar pulsations (sim- 
ulated by variable inner boundary conditions; cf. Sect. 2.3) 
and the heating of the gas caused by collisions with drift- 
ing dust grains (cf. Sect. 2.1). Concerning the interaction be- 
tween the gas and the radiation field, we have studied mod- 
els based on two different representations of the gas opacity 
adopted in the literature. On the one hand we have calculated 
models assuming a constant value of the mass absorption co- 



efficient (constant-opacity models); an approach widely used 
in time-dependent wind models. We have also, on the other 
hand, calculated models based on Planck mean absorption co- 
efficients, based on molecular data (Planck mean models). This 
approach results in more realistic pressure-temperature condi- 
tions in the photospheric layers, and lower densities in general 
(cf. Sect. 2.2). 

We have selected the sets of model parameters from a sam- 
ple where we expected that drift could be important to the 
wind structure. In the selection process we used stellar param- 
eter relations derived from observations and evolution models 
to constrain the number of independent parameters - thereby 
more clearly separating the effects of drift. Spatially our mod- 
els cover a region ranging from below the stellar photosphere 
out to about 25 stellar radii, enclosing the zones of dust forma- 
tion. 

We have compared drift models with the corresponding 
position coupled (PC; non-drift) models, and have found that 
the general changes in the spatial structure and variability in 
Paper I, caused by drift, are confirmed in the larger set of mod- 
els presented in this article. When drift is allowed dust tends to 
accumulate to the dense regions behind shocks. The gas-dust 
collisional interaction in these regions is comparatively strong 
and the drift velocity conversely low, the regions between the 
shocks are correspondingly depleted of dust. The consequences 
of the spatial relocation of the dust are larger temporal varia- 
tions of dust-related quantities (at fixed locations in the wind). 
Gas-related quantities are also affected in this process, but to a 
much lower degree. 

Several PC models in this and earlier studies show a 
(multi-)periodic variability in for instance the mass loss rate 
and the terminal velocity, whereas most drift models show ir- 
regular variations. However, we have found a couple of drift/PC 
model pairs where the drift model is (quasi) periodic, and the 
PC model is not. In particular, we mention that three of these 
drift models are Planck mean models. No periodic Planck mean 
PC models have been found here. 

Considering the changes between drift models and PC 
models - in the context of temporally averaged quantities - 
we have found significant, but not dramatic, differences for 
most models. There is, for instance, a tendency towards re- 
duced mass loss rates and a larger scatter in the degree of con- 
densation in the drift models. In a few cases the inclusion of 
drift increases the mass loss rate, possibly by providing more 
efficient conditions of wind formation. This behavior is found 
in a few constant-opacity models where the corresponding PC 
model shows a low terminal velocity and a mass loss rate at the 
lower end of the range obtained here. 

The effects of drift are, however, more pronounced in the 
Planck mean models. The winds formed in these models show 
much lower gas densities, and - as a consequence of the weaker 
coupling between the gas and the dust - show stronger effects 
of drift. In correspondence with the results in Paper I we have 
found several PC models where no drift model wind forms us- 
ing the same model parameters. All Planck mean PC model 
winds without a corresponding drift model wind share a low 
variability, in particular in the degree of condensation. In ad- 
dition Planck mean drift models tend to be more sensitive to 
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the period of the stellar pulsations than constant-opacity drift 
(and PC) models are. Concluding we note that drift adds com- 
plexity to the wind structure, both in Planck mean models and 
in constant-opacity models. It is in view of the current results 
questionable if a simple recipe exists that can reproduce the 
effects of drift using only PC models. 

We have studied how the models in this article compare 
with the mass loss rate fit formula derived by Wachter et al. 
(2002). The comparison shows that the mass loss rates are sys- 
tematically lower than this fit predicts, partly due to the use 
of different parameters in the models, partly due to the fact 
that the fit was derived for constant-opacity non-drift models. 
Planck mean PC models give mass loss rates which are lower 
than the values of the corresponding constant-opacity models 
(see above). The fit values are about ten times larger than the 
computed rates of Planck mean models, with or without drift. 
The fit formula is in its current form inapplicable to the models 
presented here. 

Although the models in this article are more advanced than 
previous models they are not yet directly comparable with ob- 
servations. One essential component, in this context, that is not 
yet included is frequency dependent radiative transfer (see e.g. 
Hofner 1999; Hofner et al. 2003). We expect that the conditions 
in the wind formation zone of more realistic (non-gray) models 
will change to somewhere in between the two cases represented 
by the two different gas opacities studied here, probably closer 
to the Planck mean models with their lower densities. 

Also worth stressing is that we so far only study the short- 
term evolution in our models (i.e. covering up to a few hundred 
pulsation periods). A study of long-term evolution, important to 
stellar evolution for instance, requires longer calculated time- 
intervals (covering thousands of pulsation periods). We intend 
to study these aspects closer in future papers. The last point 
concerns the direct microscopic influence of drift on the dust 
formation process (nucleation and growth rates). We have so far 
assumed that the stellar winds will not be significantly aff'ected 
by this effect (cf. Sect. 2.1). However, we will - for consistency 
- in the near future include these modifications in our models 
and study the dust formation, affected by drift, in more detail. 
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